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The paper that is commented by Touchette contains a computational study which opens the door 
to a desirable generalization of the standard large deviation theory (applicable to a set of N nearly 
independent random variables) to systems belonging to a special, though ubiquitous, class of strong 
correlations. It focuses on three inter-related aspects, namely (i) we exhibit strong numerical indi- 
cations which suggest that the standard exponential probability law is asymptotically replaced by 
a power-law as its dominant term for large TV; (ii) the subdominant term appears to be consistent 
with the q-exponential behavior typical of systems following (/-statistics, thus reinforcing the ther- 
modynamically extensive entropic nature of the exponent of the q-exponential, basically N times the 
(/-generalized rate function; (iii) the class of strong correlations that we have focused on corresponds 
to attractors in the sense of the Central Limit Theorem which are Q-Gaussian distributions (in 
principle 1 < Q < 3), which relevantly differ from (symmetric) Levy distributions, with the unique 
exception of Cauchy-Lorentz distributions (which correspond to Q = 2), where they coincide, as well 
known. In his Comment, Touchette has agreeably discussed point (i), but, unfortunately, points 
(ii) and (iii) have, as we detail here, visibly escaped to his analysis. Consequently, his conclusion 
claiming the absence of special connection with (/-exponentials is unjustified. 

PACS numbers: 02.50.-r,05.20.-y,05.40.-a,65.40.gd 

Before addressing in detail the Comment by Touchette [lj on our paper let us describe the physical scenario 
within which we have undertaken a possible generalization of the standard large deviation theory (LDT). A standard 
many-body Hamiltonian system in thermal equilibrium with a ther- mostat at temperature T is described by the 
Boltzmann-Gibbs (BG) weight, proportional to e"^™ = e~^ UN / N]>N , where TLn is the JV-particle Hamiltonian, and 
ft = 1/ksT. For standard Hamiltonian systems (typically involving short-range interactions and an ergodic behavior), 
the total energy is extensive. Consequently, the quantity [Hn/N] scales with N, analogously to a (thermodynamically) 
intensive variable. This is to be compared with the LDT probability P(N) ~ e~ TlN , where the rate function ri (the 
meaning of the subindex 1 will soon become clear) is related to a BG entropic quantity per particle, and plays a role 
analogous to (3[Hj\r/J\f] (we remind that, for such standard systems, /? is an intensive variable). 

If now we focus on say a <i-dimensional classical system involving two- body interactions whose potential asymptot- 
ically decays at long distance r like — A/r a (A > 0; a > 0), the canonical BG partition function converges whenever 
the potential is integrable, i.e. for a/d > 1 (short-range interactions), and diverges whenever it is nonintegrable, i.e. 
for < a/d < 1 (long-range interactions). The use of the BG weight becomes unjustified ["illusory" in Gibbs words 
Q for say Newtonian gravitation, which in the present notation corresponds to (a, d) — (1,3), hence a/d = 1/3] 
in the later case because of the divergence of the BG partition function. We might therefore expect the emer- 
gence of some function /(Hat) different from the exponential one, in order to describe some specific stationary (or 
quasi-stationary) states differing from thermal equilibrium. The Hamiltonian T-Ln generically scales like NN with 
N = N d 1 = m a/d N (with the q-logarithmic function defined as \n q z = z 1 l q 1 ', z > 0; lni z = lnz). Notice 
that (N -)• oo) N ~ N 1 '^/^ - a/d) for 1 < a/d < 1, N ~ IniV for a/d = 1, and N - l/(a/d - 1) for a/d > 1. 
The particular case a = yields N ~ N, thus recovering the usual prefactor of Mean Field theories. The quantity 
/3H N can be rewritten as [(0)Hm/(NN)]N = \j?H N /(NN)]N, where (3 = @N = l/k B f = N/k B T plays the role 
of an intensive variable. The correctness of all these scalings has been profusely verified in various kinds of thermal 
0, diffusive @ and geometrical (percolation) systems (see also 0, We see that, not only for the usual case 
of short-range interactions but also for long-range ones, [/3TLn /(NN)] plays a role analogous to an intensive variable. 
The q-exponential function e z q = [1 + (1 — q)z] T ^ (ef = e z ) (and its associated q-Gaussian [9() has already emerged, in 
a considerable amount of nonextensive and similar systems (see [Tol - [3l| among others) , as the appropriate generaliza- 
tion of the exponential one (and its associated Gaussian). Therefore, it appears as rather natural to conjecture that, 
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in some sense that remains to be precisely denned, the LDT expression e TlN becomes generalized into something 
close to e q r " (q £ TV), where the generalized function rate r q should be some generalized entropic quantity per par- 
ticle. Let us stress a crucial point: we are not proposing for long-range interactions, and other nonstandard systems, 

— r 

something like e q q with 7 ^1, but we are expecting instead 7 = 1, i.e., the extensivity of the total q-generalized 
entropic form to still hold (32|, in order to be consistent with many other related results (e.g., [si l33l - [3a |). We shall 
soon sec that this important assumption indeed appears to be verified in the model, characterized by (Q,j,S) , that 
we numerically studied in Q. 

Let us start by exhibiting that its N — >• oo LDT asymptotic behavior numerically satisfies 
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7 o_?^ > 0- This implies the existence of a generically positive finite B(x) such that 
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P(N;n/N < x)N^ 



B(x) 



N = C(x), 



(2) 



C(x) being a generically positive finite number for all values of x different from 1/2. This is indeed verified, as 
exhibited in Fig. [T] and [21 More precisely, we verify for fixed (Q,7,<5) that B(x) is unique for any given x, whereas 
C(x) is in fact a set of values, noted {Cj(x)}, with j = 1,2, . . . ,j max (the value of j ma x depends on x; for example, 
we can see that, for the illustration exhibited in Fig. (5J j max = 10 for x = 0.1). Let us emphasize that the 1/N 
correction to the power law 1 /N v in (JT]) is consistent with the total entropy of the system always being extensive in 
the thermodynamical sense. 

Let us next check the conjecture made in [2], namely that P(N; n/N < x) is, for q > 1, well approached by 
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By identifying this expansion with Eq. ([1]) we obtain 

a J (x)=B(x)[(q-l)r q ^(x)}^ 
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Since -B(a;) and {Cj(x)} are numerically known, we can easily calculate {a,j(x)} and {fg (x)} by using Eqs. Q 



and ([S]). Knowing these, we calculate a, (a 



- r u;(a;)JV 



(j = 1, 2, . . . ,j m ax) an d compare with our numerical data. We 



then bound our numerical results from both below and above (see Figs. [3] and 0] for illustrations). More precisely, 
for each value of x, we have adopted two values, noted Ci ower bound(x) and C upper bound{x), such that g-exponential 
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B=4. 57661 
B=4.576608 
B=4.57660795 
B=4.57660794 
B=4.57660793 
B=4.5766079 
B=4. 5766078 
B=4. 5766075 
B=4. 576607 
B=4. 576606 



FIG. 1: Detailed numerical verification of the conjecture given by Eq. ([2]) in one of the 10 "lines" (the bottom one, to be more 
precise) observed in Fig. [2] for x — 0.1 . This procedure enables a high precision numerical determination of B(x) for any 
chosen value of x. For a given (Q, 7, 8) model, the value of B(x) is one and the same for all the "lines" associated with a given 
value of x. Not so for C(x): indeed, for fixed x, we observe the existence of a set of values for C(x) which we note {Cj(x)}, 
with j = 1, 2, . . . ,jmax (in the present illustration j ma x = 10). The finiteness of the set {Cj(x)} here and in Fig. [2]means that 
the corrections to the N~ v power in Eq. {TJ are of the 1/N order. The finiteness that we observe (here and in Fig. [2]) in the 
slopes at the origin means that the next corrections are of the 1/N 2 order. In this example, we have run N up to 11 x 10 6 . 



upper and lower bounds for the entire set of numerical values for P(N\n/N < x) are obtained. These Ci OW er boundix) 
and Cupper bound(x) values turn out to be comparable to the corresponding set {Cj(x)} (see Fig. [21). In other words, 
we obtain the values of Gi ower bound{x), C U p per bound(x) , r q upper bound ) ^ anc j r ^° wer bound ) consistent with Eq. ([5]). 
By introducing these values in Eq. Q we obtain aq upper bound ) ^ anc j aq lowerbound ^ [x), and verify that, in the model 
studied in gj], Vx, WN, 

B(x)[(q - l)A lowerb ° und \x)]^e~ r «°^ rb ° U " d){x)N 
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We may summarize the above considerations by conjecturing that, for all strongly correlated systems which 
have Q-Gaussians (Q > 1) as attractors in the sense of the central limit theorem (see [33[), a model-dependent 
set [q > l,B{x) > 0,r q lowerbound) (x) > 0,r q upperbound) (x) > 0] might exist such that P(N;n/N < x) 
generically satisfies inequalities ([6]). In our present example, this set depends on (Q,7, S). Typical values of 
[r q lower bound) (x), r[ upper bound) (x)} are illustrated in Fig. [5] and compared with q-generalized entropic quantities. 

Touchette mentions Kaniadakis' K-logarithm and /t-exponential [361 ] as an alternative to the g-exponential and 
q-logarithm herein conjectured. Let us address this point through the definition 



ln K . z 



^ = i(ln 9 z-ln 9 i) (q = l + n) (7) 
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FIG. 2: Numerical verification of the conjecture given by Eq. ([2]), N running up to 11 x 10 6 . For fixed x, B{x) is unique, whereas 
C(x) corresponds to a set of values {Cj} (j — 1,2,..., j m ax), where jmax depends on x. We also see that the next correction is of 
the type 1/N 2 . The upper and lower values C m in and C max that are indicated by the arrows precisely correspond to the upper 
and lower bounds r ( lower bound ) ^ an( j r (" wer bound) ^ sucn that all present numerical results are within two g-exponentials, 
as indicated in Figs. [3] and 2] 
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FIG. 3: Comparison of our numerical data (dots) with a(x)e q Tq , where (a(x),r q (x)) have been calculated from (B(x),C(x)) 
by using Eqs. @ and (JSJ). The values for C{x) that have been used are those indicated by arrows in Fig. [2] Two values for x, 
namely x = 0.10 and x = 0.35, have been illustrated here. 
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FIG. 4: The same data of Fig. [3] in (g-log)-linear representation. Let us stress that the unique asymptotically-power-law 
function which provides straight lines at all scales of a (q-log)-linear representation is the (/-exponential function. The inset 
shows the results corresponding to N up to 50. 



(Notice a misprint in the definition of the K- logarithm appearing in Touchette's Comment). It straightforwardly 
follows the asymptotic series 
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The dominant term is a power-law, and at this approximation it is trivially as admissible as virtually any other 
power-law. However, we verify a highly meaningful discrepancy with the g-exponential function, namely that its 
subdominant correction is in 1/N 2 , instead of 1/N. This fact excludes the ^-exponential function as an adequate one 
for the present purpose; indeed, it cannot satisfactorily reproduce the results exhibited in Figs. [3]and|3]. 

A point remains to be discussed. The Levy-Gnedenko theorem concerns sums of infinitely many independent (or 
nearly independent, in a specific sense) random variables, whereas the 2008 Q-central limit theorem |33j concerns 
sums of infinitely many strongly correlated variables within a specific class. The first case corresponds to divergent 
standard variance, whereas the second one concerns finite Q-variance (Q = 2Q — 1; see details in [3j]). The attractors 
for the former case are Levy distributions, whereas those for the latter are Q-Gaussians. Both classes have long tails. 
For the Levy distributions, the decay is slower than l/|x| 3 and faster than l/|x|; for the q-Gaussians, any power-law 
faster than l/|x| is admissible. It is known that they have this and other relevant differences. They always differ 
excepting for an unique case, which happens to be precisely the case focused on by Touchcttc, i.e. Q = 2, namely 
the Cauchy-Lorentz distribution (named after Cauchy by mathematicians, and after Lorentz by physicists). They can 
be simply thought as having ri(x) = 0, which, as acknowledged by Touchette, is not particularly enlightening. But 
they can be also thought in a much more interesting way, namely as having r%{x) different from zero, which neatly 
illustrates the usefulness of the approach adopted in Q. In fact, it is well known that Q — 2 is a highly peculiar 
case within the interval 1 < Q < 3. For example, the anomalous diffusion coefficient in the nonlinear Fokker-Planck 
equation known as the Porous Medium Equation and discussed in [18J changes its sign precisely at Q = 2 (see also 



6 



1.00 



0.75 



0.50 



0.25 




FIG. 5: ^-dependences of the lower and upper bounds for the rate function r q (x) obtained from g-exponential fittings of the 
numerical data (see Figs. 0and[3]). The analytical curves 73/2(2) and I 5 / 3 (x) are included for comparison. The inset exhibits 
the quasi-parabolic behavior at both sides of x = 1/2. 



[HI). The fact that, for Q = 2, r± = whereas r q ^ is totally analogous to a variety of dissipative one-dimensional 
maps whose Lyapunov exponent vanishes at the edge of chaos. In such cases, the use of the nonadditive entropy S q 
instead of the BG one makes the discussion much richer since it enables a simple quantitative characterization of the 
nonlinear dynamical behavior (by generalizing the standard exponential sensitivity to the initial conditions when the 
maximal Lyapunov exponent is strictly positive to the g-exponcntial form at the edge of chaos, when the maximal 
Lyapunov exponent vanishes). This has been verified both analytically and numerically in very many cases 12|. 

Let us conclude by saying that point (i) of the present Abstract is agreeably discussed in Touchette's Comment, 
but a neat analysis of the important points (ii) and (iii) is notoriously absent in his paper. In other words, the 
g-exponential ansatz proposed in [2| for (asymptotically) generalizing the standard LDT remains (either exactly or 
approximatively: see the quantity ([2]), expected to be finite, and the inequalities |6])) as a very strong candidate for 
a wide class of systems whose elements are strongly correlated. This fact may be seen as a strong indication that, 
consistently with other results available in the literature (see 0-11 Hi EH HI] ) 5 the total entropy remains extensive (i.e., 
thermodynamically admissible) even in nonstandard cases where the BG entropy fails to be extensive. Any analytical 
results along these or similar lines would obviously be highly interesting and welcome. 
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